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Abstract: Magnetic and Rayleigh dark matter are models describing weak interactions of 
dark matter with electromagnetism through non-renormalizable operators of dimensions 5 and 
7, respectively. Such operators motivate the existence of heavier states that couple to dark 
matter and are also charged under the electroweak interactions. The recent hints of a gamma- 
ray line in the Fermi data suggest that these states may be light enough to be produced at the 
LHC. We categorize such states according to their charges and decay modes, and we examine 
the corresponding LHC phenomenology. We emphasize unconstrained models that can be 
discovered in targeted searches at the upgraded LHC run, while also enumerating models 
excluded by current data. Generally, models with SU^(2)-singlet states or models where the 
charged states decay predominantly to tau leptons and/or gauge bosons are still viable. We 
propose searches to constrain such models and, in particular, find superior performance over 
existing proposals for multi-tau analyses. Finally, we note several scenarios, especially those 
dominated by tau final states, that cannot be probed even with 300/fb at LHC 14, motivating 
the further refinement of tau lepton searches to improve sensitivity to such final states. 
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Figure 1: Estimates of model exclusions and discovery reach for fermions (red/upper bar) and scalars 
(blue/lower bar). Solid color indicates an excluded model, while light shading means that it can 
be discovered at 5cr with < 300 fb~^ at LHC14. The top (bottom) plot shows models with flavor- 
dependent (-independent) couplings. Solid black indicates unstable dark matter. Models are described 
in Section 2 and constraints/discovery reach in Sections 3 and 4. 
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1. Introduction and summary 

The weakness of interactions is often understood in field theory as a sign that the corre- 
sponding operators are irrelevant. Consequently, the "darkness" of dark matter (DM) may 
be naturally interpreted as a consequence of DM having only irrelevant interactions with 
light, and more generally with the electroweak gauge bosons. The lowest- dimension irrele- 
vant operator is the dipole operator 

/:dipole = (^) XV^.S'^'^X + h.C. (1.1) 

This operator requires the existence of two separate states x %*, or it otherwise vanishes. 
In its absence (as is the case for a Majorana particle), the next most-relevant operator is 
the Rayleigh operator, which can couple to both hypercharge and SU^(2) gauge bosons. We 
parameterize these interactions as 

/:Rayieigh = ^ {xx (cos O^B^.B^'^ + sin e^l^W^uW^''') (1.2) 

+ i X75X (cos 0^ B^^B^" + sin 6^ TtW^.W^'") } . 

Both dipole and Rayleigh operators give rise to the scattering, annihilation, and production 
of DM. The impact of the magnetic operator in a variety of observations has been the subject 
of several past works [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11]. The Rayleigh operator was studied in 
more recent dedicated analyses [12, 13, 14, 15], although older studies already appreciated its 
relevance for low energy experiments [2]. 

In the past year, a new impetus to these considerations was given by the observation 
of an excess of gamma-ray events from the center of the galaxy in the Fermi satellite data; 
see refs. [16, 17, 18, 19, 20] for early investigations, refs. [21, 22, 23, 24, 25, 26] for more 
detailed follow-ups, and refs. [27, 28, 29, 30, 31] for constraints on continuum emissions. 
Given the non-renormalizable nature of the above interactions and their increasing strength 
at high energies, the LHC is a promising place to look for their signatures. This was already 
considered to some extent in [15], with a general operator analysis, but it proves difficult to 
say anything definitive in the most general case. More recently, in ref. [32], a subset of the 
authors considered a simple renormalizable theory which gives rise to the interactions (1.1) 
and (1.2). Aside from the DM state, this theory includes a charged scalar {(p) and a charged 
fermion with a Yukawa interaction of the general form 

CDXi^X'P- (1.3) 

In ref. [15] it was shown that, if the theory is perturbative, the charged states should have 
masses below approximately 300 GeV. Since they are charged under the electroweak group, 
pair production of these states occurs at the LHC. The purpose of this paper is to consider 
the phenomenology of such electroweak production, and to highlight the LHC searches in 
which these new states may be discovered. 
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We emphasize that, while the current work was motivated by the specific model of [32], it 
is relevant more generally to the production of 0(100 GeV) weakly charged states. As such, 
the bounds presented and searches proposed in this paper are applicable more generally to 
many other models. Indeed, several other models were proposed in connection with the Fermi 
line [33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49], many of which feature 
new weakly interacting states at the electroweak scale whether explicitly or implicitly. Other 
collider studies related to the Fermi line can be found in refs. [50, 51], but our work studies 
a wider range of charge assignments and final states, assessing their status and observability. 

In this paper, we consider models where the charged fermion, ij;^ and the charged boson, (/p, 
carry only electroweak charges and cannot be strongly produced. If and (p carry color charge, 
then loop-induced coupling of dark matter to gluons results in a bound of m^p ~ > 500 
GeV from LHC monojet searches and the XENONIOO direct detection experiment [52], in 
contradiction with the requirement from perturbativity that ~ m^p < 300 GeV. 

Our findings are presented in Fig. 1, which sums up the current experimental constraints 
on each model considered and the discovery prospects for LHC14. We find that generic 
models with new states decaying into gauge bosons and light-flavor leptons (e//x) are almost 
entirely ruled out up to 300 GeV by current multilepton analyses. Models with multiple r 
leptons in the final state, and/or with SU^(2) singlet charge, are less constrained but are 
mostly within reach of LHC14. Throughout the paper, we use a benchmark of 300 fb~^ to 
determine whether a model can be discovered, although in many instances, considerably less 
luminosity is required for complete discovery in the 100-300 GeV window. There are, however, 
a few models with very challenging signatures (such as t^t~ + final states, disappearing 
charged track signatures [53, 54, 55, 56, 57, 58, 59], or entirely hadronic decays) that may not 
be accessible at LHC 14, even with large integrated luminosity. 

Table 1 outlines the classes of searches most relevant for the exclusion and discovery 
of the different models we consider. The most constraining searches to-date at the LHC 
are multilepton-f^Er for final states with e and /i, and same-sign dilepton searches with 
an additional hadronic tau for multi-tau final states. These analyses will continue to be 
important at 14 TeV for discovering or constraining the models we study, although we suggest 
in some instances how these searches can be optimized for particular final states by changing 
kinematic cuts, requiring additional hadronic tau tags, and using resonance reconstruction. 
In particular, we note that same-sign dilepton + tau searches are ideal for discovering tau-rich 
final states, and for certain benchmark models (such as a Type IV Two Higgs Doublet model), 
the searches we propose allow discovery at much lower luminosities than other proposed 
analyses [60], and at fixed luminosity they improve the signal significance by factors oi ^ A 
in 3t final states and ^ 2 in At final states (details can be found in Sections 3.3 and 4.3). 

The paper is organized as follows: readers who are less interested in the details can 
consult Fig. 1 and the conclusions in Section 5 for a summary of the results of this paper. 
Due to the large number of models, we focus throughout the paper on the most constraining 
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Collider search 


Model 




Sterile neutrino model {^p coupled to e, /i) 




Odd scalar model (e and ji in final state) 
Odd fermion model (Type I 2HDM, triplet) 


Same-sign dilepton + hadronic r 


Odd scalar model (doublet, r in final state) 
Odd fermion model (Type IV 2HDM) 


Same-sign dilepton + dijet resonance 


Odd fermion model (triplet) 


Disappearing charged tracks 


Stable model 


No distinctive signature 


Sterile neutrino model {if coupled to r) 
Odd scalar model (singlet, r in final state) 



Table 1: Overview of collider searches relevant for constraining the models in Section 2, whether with 
current or future data. Details of the searches are given in Sections 3 (current bounds) and 4 (future 
searches) . 



signal(s) for each scenario, and we outline the classes of searches relevant for the exclusion 
and discovery of each in Table 1. We begin the paper by categorizing the different models in 
Section 2. In Section 3, we discuss the bounds on such models from current LHC searches, 
and outline proposals and prospects for discovery at 14 TeV in Section 4. Details related to 
our analyses can be found in the Appendices. 

2. Models for the new charged states 

Our motivation is the study of charged states coupled to DM that generate DM-photon 
couplings as in Eq. (1.3). The interaction in Eq. (1.3) does not, however, completely fix 
the phenomenology of the model, since other interactions may allow the states Lp and 
to decay. There are several phenomenological classes of models, and we discuss the four 
different possibilities in this section. The implications of each of these scenarios for LHC 
phenomenology are addressed in the subsequent sections. 

The most straightforward and minimal resolution to the problem of new stable charged 
states is already present in the model itself. The interaction Eq. (1.3) respects a ^7(1)^ 
symmetry under which the DM and messengers carry charges; this symmetry forbids any 
other coupling of ijj by themselves to any current made of SM fields, and implies the 
existence of at least one other stable particle in addition to the dominant dark matter x- For 
an SU^(2) multiplet, the electromagnetic corrections make the electromagnetically charged 
components of the multiplet slightly heavier. This is similar to the situation occurring in 
models of anomaly- mediated supersymmetry (SUSY) breaking [61, 62], in which the wino 
state is the lightest supersymmetric particle (LSP). The charged components have a long 
lifetime and decay to the neutral components through the emission of an off-shell W^-boson. 
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Table 2: Possible electroweak charge assignments for ip/cp are hsted in the first column. Higher mul- 
tiplets and larger hyper charges are also allowed, but are mostly excluded due to their high production 
cross sections. In the second column, we list the allowed renormalizable interactions with SM fields 
when the scalar cp is Z2 odd {cr^ are the SU^(2) generators). We ignore operators which are equivalent 
up to a field redefinition (such as (fx^^ for charge (1,-1)). In the third column we list the correspond- 
ing operators when the fermion ip is Z2 odd. We use V{(f,h) to indicate a generic two-scalar potential 
coupling (f to the SM Higgs field. 

This possibility is particularly intriguing since it results in more than one type of stable 
neutral particle. Most of the DM is typically constituted by the x state since it is lightest, 
but some relic abundance of the neutral components of ^ and may be present. We discuss 
this model and the cosmological implications in Section 2.1. 

Another possibility, with markedly diflFerent phenomenology, is to endow one of the mes- 
sengers {ip or ^) with some decay mode into SM particles. However, at least one messenger 
state together with x must be odd under a Z2 symmetry; otherwise the DM candidate, is 
itself rendered unstable due to the interaction Eq. (1.3). We therefore draw the distinction 
between models according to which of the charged particles, scalar or fermion, is Z2 odd (i.e. 
has no direct decay into only SM particles). 

Another discrete choice that must be made is the electroweak charge carried by the 
messenger fields. The representations we consider in this paper are shown in Table 2. Each 
charge assignment results in diflFerent renormalizable couplings between ip/^ and SM fields-"^. 
We now present each of these possibilities in turn. 

2.1 Stable model 

The first model we consider is one where the lightest neutral component of cp and/or ijj is 
stable (in addition to the dark matter %). For concreteness, we consider ip^ as the lightest 
component, although our results extend easily to the scalar case. The neutral fermion ip^ is a 
thermal relic with abundance determined by its annihilation to XX through the Yukawa cou- 

^ Larger SU^(2) multiplets and hypercharges are possible, but they tend to produce similar final states and 
the production cross sections are greatly enhanced over those in Table 2, making constraints much more severe 
and discovery more straightforward. 
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pling in Eq. (1.3)^. The interaction strength, A, is required to be large in order to account for 
the observed Fermi gamma-ray hne, and this ensures that does not overdose the universe. 
Nevertheless, we wish to determine whether the ip^ relic is still large enough to violate exper- 
imental constraints on DM, particularly direct detection bounds that strongly constrain the 
DM fraction coupling to the Z. We start by determining the t/j^ relic abundance. Assuming 
5- wave- dominated annihilation, the thermally-averaged tjj^ annihilation cross section is 

|A|4 2m| - ml 




We calculate the relic abundance Jl^ following the derivation of [63] , and subsequently 
compute the relic number density of relative to that of a weakly interacting massive 
particle (WIMP) of the same mass which saturates the DM component. We find that the 
number density of tp^ is suppressed relative to DM of the same mass by a factor of 5 x 10~^ 
to 10~^, depending on the precise values of A, m^, and m^^. 

For m^, rrii^ ^ 130 — 300 GeV, which is our range of interest, a thermal relic with a 
(2, —1/2) electroweak charge is constrained by XENONIOO to have a density that is < 10~^ 
times that of dark matter [64]. As a result, essentially all of the models with stable ip^ 
are ruled out if they couple to the Z, even though they make up only a tiny fraction of 
DM. Therefore, the most viable models with stable ip^ are those where the messenger fields 
have zero hypercharge, and as a consequence, the neutral states do not couple to the Z. 
The authors of [51] found a similar result, even though annihilation proceeded through a 
different channel in that case.^ 

When tp^ is stable, the charged components in the SU^(2) multiplet decay into it. The 
mass splitting between charged and neutral particles in the multiplet is generated radiatively 
after electroweak symmetry breaking, and for ^ 100 GeV, the splitting is 0(100 MeV). 
For a splitting induced entirely by electroweak gauge boson loops. Am = — m^- > m^r, 
and production of ip^ produces V^^, along with a soft pion. Tree-level modifications to the 
splitting can increase the lifetime and can give Am < m^^, leading to a decay to a soft lepton 
and neutrino. We show the decays in Fig. 2. 



2.2 Odd scalar model 



When the charged scalar messenger is odd under the Z2, the fermion is Z2 even and it can 
mix with the leptons of the SM, rendering it unstable. The relevant operators are shown in 
the second column of Table 2 for the different charge assignments. The mixing is either a 

^The relic densities of cp and ip are also affected to a small degree by the messengers' interactions with 
electroweak gauge bosons, akin to the Higgsino in the MSSM. 

^If the Dirac fermion is split into a pseudodirac state, or the complex scalar is split into non-degenerate real 
and imaginary pieces, the Z-coupling can be suppressed [65, 66, 67]. However, this requires a higher dimension 
operator, implying additional fields. For economy, we neglect this case. 
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Figure 2: Decays of the charged state into the stable neutral state ip'^. The left (right) decay 
dominates when the mass difference — m^o > (< m^). 

direct bilinear mixing or it is achieved after the Higgs boson develops a vacuum expectation 
value (VEV)^. We discuss in some detail the (2, -1/2) charge assignment, and comment on 
the other possibilities below. 

The possible interactions are 

jC D Mf^^^ + A^£/ie^ + yt/;he^ + m^^^ (2.2) 

where h is the SM Higgs, £ (e^) is any one of the left (right)-handed leptons, and flavor indices 
are suppressed. After electroweak symmetry breaking, this results in a mixing between the 
SM charged leptons and the heavy SU^(2) doublet fields t(j = ('0^,'^/^") and ij;^ = ('^/^^^, '^/^+) , 

jC D Mf^V^° + Mf^-^+ + yv ^-e^ + ee^ + mu^^^ + me^+. (2.3) 

The neutral component is easily unmixed by defining 

= ^ (Mf^° + mv) . (2.4) 

This introduces a coupling between the WIMP (x) and the scalar messenger [ip) to the light 
neutrinos v through the Yukawa vertex Eq. (1.3). It allows the neutral scalar component 
to decay via ^ X + ^- The charged fermion component similarly mixes with the SM 
electrons and allows for the decay of the charged scalar component if'^ X +6^- The mass 
mixing of the fermions with the SM leptons allows the fermions to decay through W^, Z, or 
h emission into SM leptons. These decays are shown in Fig. 3. We note that the mixing must 
be fairly small; otherwise, it introduces a new annihilation channel of DM into leptons that 
can easily dominate over the annihilation into photons [32] . Similarly, couplings to electrons 
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Figure 3: Decays of the new charged states when the heavy fermions mix with SM leptons. On the 
left, the scalars decay into the WIMP candidate, Xi a SM lepton. On the right the fermions decay 
directly into a gauge-boson and a SM lepton. 

and muons have to satisfy y < 10~^ to avoid excessively large contributions to anomalous 
magnetic dipole moments. 

If the charged fermion is a singlet under SU^(2), the mixing can be done in an analogous 
fashion, although mixing now occurs with the right-handed electrons. The main difference is 
that there are no heavy neutral states in this case. Finally, in our discussion so far we have 
neglected the issue of flavor. Flavor-changing neutral current constraints place very strong 
bounds on models where ^ mixes with different flavors of leptons. We take a phenomenological 
approach and assume that couples exclusively to one flavor of lepton (or, alternatively, that 
the couplings are aligned with the SM Yukawa matrices [68, 69, 70]). If a signal is seen in 
any channel, it would be important to determine how the result is compatible with flavor 
observables. 

2.3 Odd fermion model 

When '0 is odd under the Z2 symmetry, operators exist that allow to decay into SM 
fields. Most such models couple to the SM through the scalar potential including the Higgs 
boson, V{h^Lp). Depending on the terms in the potential and the quantum numbers for (/p, 
components of Lp can mix with the Higgs doublet components /i, and the neutral component 
can develop a VEV. A VEV for Lp is constrained since it induces a Dirac mass term mixing x 
with '0 through the Yukawa coupling X{(f)xi^' If ^ violates custodial SU(2) symmetry, more 
stringent constraints on {(p) arise from bounds on the p parameter [71, 72]. In the following 
sections, we consider models where (p mixes with the SM Higgs as well as models where it 
couples directly to fermions. 



^ After electroweak symmetry breaking, the two are equivalent up to field and coupling redefinitions. 
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Lp mixes with the Higgs boson 



The scalar Lp can develop a VEV in one of two ways: through a negative mass-squared 
term or through a linear term. The former is challenging to realize in a model, since the 
VEV and mass of are connected by the quartic coupling in the Lp potential, X^. The precise 
relation is = \^ (9^)^/2 <^ even if A^^ is large. Doublets below mz/2 are generally 
excluded, so this is not a valid way to induce a ip VEV. The other possibility is to have a 
source for and we now consider several models with this feature. 



(2,-1/2): LP and h comprise a Two Higgs Doublet Model (2HDM). Terms of the form 
(as well as induce a source and mix the scalars. If the mixing term is of the form 

Mmix^*^5 the induced VEV for (p is {(p) ^ (/i^-^//i^) t', where /i^ is the coefficient of the 
quadratic term for ip. Diagonalizing the masses gives the mixing angle sin a: 

sma^^^. (2.5) 

In the mass basis, there are three heavy scalar states: using conventional 2IIDM terminology, 
these are the heavy CP-even and CP-odd scalars, H and A, and the charged scalar H^. 

When {(p) ^ 0, the lightest DM mass eigenstate acquires a doublet fraction / ^ X{(p) / \/2{m^- 
m^) through the Yukawa coupling Eq. (1.3). The doublet fraction of DM is constrained to be 
/ < 0.2 — 0.45; otherwise, XX ~^ W"^W~ would violate constraints from the Fermi continuum 
gamma ray spectrum^ [27, 28, 29, 30, 31]. With = 130 GeV, at the weak scale, and 
the Yukawa coupling A being sufficiently large to generate appropriate Magnetic and Rayleigh 
operators, this implies that 

V Xv /0.4\ /3\ /100GeV\ 
tan^ = — = — > 14 — - , 2.6 

(W / i^^P - ^x) \ f J - ^xJ 



where v = 246 GeV is the SM Higgs VEV. Therefore, if (p develops a VEV, {(p) is constrained 
to be much smaller than the SM Higgs VEV. The constraints on (p)) from (2.6) are satisfied 
as long as /imix <^ j^cp- In this scenario, the mass eigenstates strongly correlate with the 
interaction eigenstates h (SM Higgs boson) and p) (heavy doublet). 

A possible problem with this model is that, considering all interactions allowed by gauge 
invariance, the heavy scalars generate flavor-changing neutral currents that are ruled out by 
many orders of magnitude. In viable 2HDMs, the scalar sector is charged under an additional 
Z2 symmetry [73], ensuring that only one particular Higgs couples to each of down quarks/up 
quarks/leptons. These are described by the different "Types" of 2HDM, and can also be 
described by models with Yukawa alignment due to a minimal-ffavor-violating hypothesis 



^The most conservative bound comes from background saturation, while a shape analysis improves con- 
straints on the annihilation cross section by an order of magnitude. 
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Figure 4: Decays of the heavy scalars (top) and fermions (bottom) in the Type I Two Higgs Doublet 
model. Scalar branching ratios depend sensitively on which final states are kinematically allowed. 



Figure 5: Decays of the heavy scalars (left) and fermions (right) in the Type IV (leptophilic) two 
Higgs doublet model. Scalar branching ratios depend sensitively on which final states are kinematically 
allowed. 

[68, 69, 70]. In what follows we focus on the most unconstrained 2HDMs: Type I, in which 
(f has no tree-level Yukawa couplings to SM fermions; and Type IV, in which (f couples 
exclusively to leptons and h couples exclusively to quarks. 

In the Type I model, the light CP-even Higgs is almost exactly SM-like due to the small 
mixing and SM- Yukawa structure. The heavy CP-even scalar H decays through mixing with 



,//^+ H 




DM 



the SM Higgs to 66, W^W~, ZZ, and hh when kinematically ahowed, while the CP-odd 
scalar A decays predominantly to hZ and bb. The charged heavy Higgs, H^^ decays to tb 
when allowed, and ru otherwise. The Z2-odd fermion can decay to a vector boson and 
the WIMP, X? due to the mixing induced by {(p). It can also decay to hx due to the mixing 
between (p and the SM Higgs. At low masses, the decays of to DM through off-shell gauge 
bosons dominates, while for > 220 GeV, ip decays predominantly to on-shell bosons. The 
decays are shown in Fig. 4. 

For Type IV models, the scalar field (f couples exclusively to leptons, which obtain 
their mass through {(p)^. Consequently, the r Yukawa coupling is 0{1). The light Higgs is 
SM-like, although the coupling to r is modified. Somewhat counterintuitively, the h t~^t~ 
rate can actually be enhanced by an factor because the large Yukawa coupling of r to 99 
compensates for the small mixing between h and (f. The decay modes of the heavy Higgses 
i7, A, and are almost exclusively to r final states. The Z2-odd fermion ip decays via 
ip^ r^UrX 2tnd ip^ t^t~x- The decays are shown in Fig. 5. 

We briefly comment on 2HDM of Types II and III, where (f couples predominantly to 
down-type quarks. These models predict enhanced couplings of heavy Higgses to b quarks 
at large tan/3. Generally, heavy Higgses are excluded by 6 ^ 57 constraints in the range of 
interest of our models (m^^i < 300 GeV) for all values of tan/3 [75]. Therefore, we do not 
consider these models further. 

(3,0): The term h^ip^a^h appears in the potential (cr^ are the SU^(2) generators), mixing 
the neutral component of Lp and h. Unlike a 2HDM, however, there is no direct coupling of ip 
to the Z, leading to a tree-level correction to the p parameter when {p) ^ 0. This means that 
the size of the mixing must be small {{p) < 8 GeV). The lack of a tree-level coupling of p to 
SM fermions means that it decays through Higgs mixing or, at high masses, to longitudinal 
W"^W~ and ZZ, and to hh. The fermions ip^/ip^ decay predominantly through the scalar 
mixing to W^/h + x- 

(3,-1): The term {eh)ip^(T^h appears in the potential. The symmetries also allow a Yuakwa 
coupling {€£)p)*^a^£ {£ is a SM lepton doublet), which would induce an unacceptably large 
Majorana neutrino mass. This possibility is therefore ruled out unless some symmetry forbids 
the p)££ interaction. As in the (3, 0) scenario, the mixing between p^ and h must be small 
((^) ^ 8 GeV) to avoid large corrections to the p parameter. The p) and tp decays in this 
model are similar to the Type I 2HDM discussed above, where the p) decays predominantly 
to SM gauge bosons and Higgs bosons when kinematically allowed, and to r and b final 
states through the Higgs mixing otherwise. Additionally, there exist doubly charged states 
p)^^ {^p^^)^ which decay to W^W^ (W^W^ + x) fi^i^l states. This is an interesting and 

^We use the convention of 2HDMC [74]; some sources invert the ordering of Types III and IV or use the 
alternative terminology of Types X and Y. 
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uncommon signature which improves the detection prospects as we discuss below in Section 4. 

if couples to fermions 

It is possible that the charge of ip does not allow it to mix with the SM Higgs. An example 
is the charge (1,-2); the only renormalizable interaction leading to Lp decay is 

/:D(^eV, (2.7) 

where are right-handed SM leptons. The scalar (^^ is a particle of charge —2 decaying into 
same-sign dileptons, giving a striking LHC signature. Other interactions are allowed with 
more exotic charges, particularly when non-renormalizable decay modes are allowed. Such 
interactions are beyond the scope of this paper, but we note that when ip is sufficiently long- 
lived and decays on detector length scales, the model becomes similar to the stable models 
we considered above in Section 2.1. 

2.4 Sterile neutrino model 

The final model we consider departs from the minimal models described above. It contains an 
additional light state with the quantum numbers of a sterile neutrino. So far we considered 
the DM state to be made up of at most two Weyl fermions, xi X2^ as in Eq. (1.1). It 
is natural to try to extend the model to three states which behave like three right-handed 
neutrinos, and so we add a third state, X3- We consider a mass structure which results in a 
pseudo-Dirac pair, xi and X2, together with a single Majorana fermion X3, 

MxiX2 + \rnixixi + 5^2X2X2 + ^rnsXsXs- (2.8) 

Introducing a heavy Dirac pair of SU^(2) doublets, ijj and ijj^ of mass Mf and a new scalar 
doublet ip with mass Mg, the Yukawa coupling is 

C D Xi^^xi + A^V^ VX2 + h.c. (2.9) 

These terms combine into a Yukawa coupling between the Dirac pairs (V^,'?/^'^) and {xi^X2)^ 
generating the operators (1.1) and (1.2). 

We consider a Z2 symmetry under which V^, V^*^, xi and X2 are odd. Then, we can write 
a Yukawa term between ip^ the normal lepton doublets, and the sterile xs? 

£ D A3£(/PX3 + h.c. (2.10) 

This allows the scalar Higgs to decay on shell to a SM lepton and the sterile state X3- It 
also allows the heavy leptonic doublets, ^ and to decay to a WIMP, a SM lepton and the 
sterile state xs through an off-shell scalar ip. These decays are shown in Fig. 6. 
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Figure 6: Decays of the heavy charged states when a hght sterile state, xa, is present. On the left the 
scalar decays into a SM lepton and the sterile state. On the right the fermion decays into the WIMP 
candidate ^ SM lepton, and the sterile state through the scalar. 

With the interactions (2.9) and (2.10), the model has an additional Z2 symmetry under 
which '0, (f, and xs odd. This symmetry can be weakly broken, for instance by mixing 
with the SM neutrinos through the operator Mxs- this case xs decay through an 
electromagnetic dipole moment transition with the SM neutrinos as is usually the case for 
sterile neutrinos [76]. If m^g > GeV then it decays before Big-Bang nucleosynthesis, and 
cosmological constraints do not apply. On the other hand, if the symmetry is unbroken 
then X3 is a stable thermal relic that is a component of dark matter, and its interactions 
are bounded by overclosure constraints. This component, xs^ annihilates to leptons through 
effective operators like A3(^X3)(^X3) V^^- For electrons ^ = e, LEP monophoton searches 
constrain Ag/m^ to be small enough that the xs I'elic abundance overdoses the universe 
unless > 20 — 50 GeV, depending on the Lorentz structure of the interaction [77]. If xs 
is a WIMP with a thermal relic abundance, it must therefore either have a weak scale mass 
or couple predominantly to muons and/or taus. 

Alternatively, xs can be extremely light or massless, in which case it contributes to the 
relativistic degrees of freedom at the times of Big-Bang nucleosynthesis and recombination. 
Its contribution depends on the temperature at which it decouples from the thermal bath, 
which depends in turn on the magnitude of A3 and the lepton species to which it couples. In 
the most constrained case, xs couples to electrons with 0(1) coupling, and it contributes a 
full effective neutrino degree of freedom (A^eff)- This is consistent with the current observed 
value [78] of A^eflf- The predicted contribution to A^gff drops significantly if A3 <C 1 or xs 
couples dominantly to taus. 
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Figure 7 : LHC production cross sections plotted against the mass of the new messenger electroweak 
state for (left) 8 TeV and (right) 14 TeV. Solid and dashed lines show cross sections for electroweak 
charges (2,-1/2) and (1,-1), respectively. 

3. LHC constraints 

Searches at the LHC for new electroweak states are hampered by the large QCD and top 
backgrounds. As a result, the only models constrained by existing LHC searches are those 
with multiple leptons in the final state, often accompanied by missing energy, New 
electroweak particles decaying to tau leptons are generally less constrained than their light- 
flavor counterparts, as are models with only gauge bosons and 0t in the final state. The 
production cross section of the new electroweak states also plays a role, as can be seen in 
Fig. 7: SU^(2) doublets or higher representations have large production cross sections and 
are more readily constrained, while SU^(2) singlets easily evade experimental bounds. 

In the following sections, we examine the collider bounds for each of the models discussed 
in Section 2. The classes of searches relevant for each model were outlined in Table 1. We find 
that, generally, the most relevant searches are the dilepton+^T, trilepton+^x, and same-sign 
dilepton analyses. For details of the simulations used in our analyses, we refer the reader to 
Appendix A. 

3.1 Stable model 

In models where the lightest neutral component of the new charged matter is stable, collider 
bounds are extremely weak. The phenomenology of such models has already been discussed 
in the literature [53, 54, 55, 56, 57, 58, 59], and there has been a recent paper applying these 
searches to a particular model of the Fermi line [51]. As a result, we sketch only a few relevant 
details here. 

The lightest neutral component of leaves missing energy in the detector and can, 
in principle, be the subject of monojet searches. However, the production cross section is 
electroweak in strength and is typically several orders of magnitude (at least 3-5) below the 
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current bounds [51]. Pair production of the charged components, (/^^/^^, can give a more 
striking signature. When the mass spHtting between charged and neutral states is ~ m^r, the 
charged particles travel some finite distance through the detector (~ 10 cm) before decaying. 
The signature consists of charged particle tracks that disappear midway through the detector, 
and monojets and monophotons can be used as a trigger [53, 54, 55, 56, 57, 58, 59]. In 
principle, there is no SM background that gives this signature. There is, however, significant 
contribution from SM combinatoric background where tracks are mis-reconstructed to give 
the illusion of a disappearing charged particle. 

ATLAS has performed a search for anomaly-mediated SUSY breaking via the direct 
production of winos, which are SU^(2) triplets displaying such a disappearing charged track 
signature. They rule out m ± < 100 GeV when the mass splitting between the charged and 
neutral state arises solely from electroweak gauge boson loops, as is expected in a minimal 
model [79]. Therefore, triplet '0 with no tree-level corrections to the mass splitting are ruled 
out up to the lower limit of our window of interest, while other charges have no constraints 
due to shorter track lengths. 

3.2 Odd scalar model 

When ^ is odd and mixes with SM leptons, the components of decay to u/i together with 
a boson W^/Z/Zi, with branching fractions determined by the SU^(2) and Uy(1) charges of 
Signal events can therefore be rich in leptons and missing energy, giving good prospects for 
discovery at the LHC. Indeed, most scenarios with tp decaying to leptons plus gauge bosons 
are already ruled out. The only models still allowed in the < 300 GeV range are those with 
small production cross sections (such as SU^(2) singlets) and those decaying predominantly 
to r leptons. Although constraints on r final states are weak, triplets and higher multiplets 
decaying to taus are ruled out by same-sign dilepton -|- hadronic tau searches due to the large 
production cross section. There are currently no constraints on the scalar decay (/9 ^ x + ^/^ 
because of the small production cross section, as well as the kinematic phase space suppression 
resulting from (f decay into dark matter, leading to low acceptances for collider cuts. 

Since the phenomenology of the mixed '0-lepton models depends on the gauge charges 
and the flavor structure, we detail each type of charge and flavor coupling in turn. We refer 
the reader to Appendix B for the fermion branching ratios of each charge. 

(2,-1/2): The doublet fermion consists of charged and neutral components, ip^. The 
neutral components decay exclusively to W^'^^^, while the charged component ip^ decays 
predominantly to Z£^^ with branching fraction ranging from 0.6 — 1 depending on m^. LHC 
production of ip^ip~ ^ ip^tp^^ and tp^tp^ leads to trilepton final states when one of the gauge 
bosons decays leptonically. ATLAS trilepton searches at 8 TeV with 13 fb~^ [80] give 

the strongest bound on the model when £ = e, /i: associated production of ip^ip^ is ruled out 
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for m^o < 350 GeV. 

If '0 decays predominantly to rV, the final state for i/j^i/;^ production is W^W~t+t~ and 
W^Zt+t", leading to same-sign dilepton + tau signatures. The strongest constraints come 
from the CMS search for same-sign dilepton plus a hadronic tau at 8 TeV with 9.2 fb~^ [81], 
which constrains m^± = m^o < 240 GeV. 

(1,-1): There exists one new charged fermion, which decays dominantly via W^u^ 
and with secondary decay Z£^. The states are pair-produced, and the most common 

LHC signature is W^W~ + which is challenging to disentangle from the large W^W~ 
background. The W^Z + £+ 0t final state is more promising, but current LHC multilepton 
analyses do not yet constrain this scenario as the fermion pair production cross section is 
small for a singlet, and most 4-lepton analyses either veto all Z bosons in an event or allow 
for 2 Zs, leading to much larger backgrounds. 

If ip decays dominantly to r, the most promising final state is WZ+r + where the 
gauge bosons decay leptonically and the r hadronically. Due to the small cross section, cur- 
rent bounds from same-sign lepton + r searches do not apply, and the prospects do not 
improve significantly at 14 TeV. 

(3,0): Analogously to the doublet model, the neutral components of the triplet decay al- 
most exclusively to W^^^, while in contrast with the doublet, charged components decay 
dominantly to W^u. As there is no Z coupling to i/j^ip , the strongest constraint comes from 
pair production of ip^ip^^ giving trilepton+^T final states when the two W bosons decay 
leptonically. The ATLAS trilepton search [80] rules out such a scenario when ^ = e, /i and 
m^± = m^o < 350 GeV. 

If decays dominantly to r, the CMS same-sign lepton plus hadronic tau search [81] 
constrains the final states i/j'^il;^ W^W^r^u. The masses m^± = m^o < 330 GeV are 
ruled out. 

(3,-1): The triplet fermion consists of doubly charged, charged and neutral components, 
^ , and ip^. Similar to fermions with charge (3,0), the fermions decay dominantly to 
W + lepton or neutrino. The pair production cross section is the highest among the possible 
electroweak charges, thus it receives the strongest constraint. For ^ = e, /i, m^± < 500 GeV 
is ruled out, while for £ = r, m^± < 400 GeV is ruled out. 

3.3 Odd fermion model 

When '0 is odd under the Z2 symmetry stabilizing dark matter, (f can decay promptly into 
SM states. The phenomenology is then dictated by the couplings of (/;, which can interact 
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either directly with SM fermions, or indirectly through mixing with the SM Higgs boson. 
When if mixes with the Higgs, its phenomenology is largely dictated by the allowed Yukawa 
couplings to SM fermions. If ip has no tree-level couplings to SM fermions, the models are 
unconstrained by current LHC data for — 100 — 300 GeV, with the exception of states 
with charge (3, —1), which are excluded for — 170 — 210 GeV. On the other hand, if Lp is 
responsible for giving mass to the leptons (Type IV 2HDM), it is constrained to have a mass 
> 180 GeV. Multilepton signatures of more general 2HDMs have been considered in [82]. 

Finally, we consider a scenario where a doubly-charged scalar (f decays directly into 
two same-sign SM leptons through a Yukawa interaction; this scenario is excluded completely 
in the window of interest unless p decays dominantly into taus, in which case — > 110 
GeV for a 100% branching ratio into taus. For — ^ 100 GeV, a sizeable branching fraction 
into taus (~ 97%) is still allowed. 

p mixes with the Higgs boson 

(2, -1/2) Type I 2HDM: After mixing with the SM Higgs, the heavy scalar eigenstates i7, 
A, and as well as the fermions '0, are directly produced at the LHC via the elect roweak 
gauge interactions. LEP has ruled out mjj± < 80 GeV in searches for ru and cs [83], and this 
remains the strongest bound to-date. Searches for direct production of H^H~ r^r~iyri^r 
are difficult due to the challenging final state (see Section 3.4 below for a similar model). The 
lack of a direct coupling between the top and H"^ renders searches in t ^ H^b ineffectual, 
while for m^± > m^, the dominant decay is to tb; the enormous top backgrounds limit the 
viability of this search. 

The heavy neutral Higgs bosons are produced in association with one another at the LHC: 
pp Z^*) AH. In the kinematic regime probed at LEP, both A and H decay to 66, and 
searches for fully hadronic decays of Higgs bosons exclude the masses up to 70 
GeV in the 46 final state [84]. This search channel suffers from large hadronic backgrounds 
at the Tevatron and LHC and there are no stronger bounds. At higher masses {mn ^ 160 
GeV), the H -^WW mode dominates, although tt background swamps the signal. The most 
promising search region is mn ^ 250 GeV, where the H —> hh and A ^ hZ modes dominate, 
but the cross sections are too small to have been probed at LHC8. 

The fermion ip also appears in the spectrum and can be pair-produced, typically decaying 
through the mixing induced between and x' ~^ + x '0^ Z + x- The branching 
ratios of '0 are given in Appendix C. When > 220 GeV, the production of the fermions 
results in diboson + signal. The leading ATLAS and CMS searches are not yet sensitive 
to i/j production in this regime. 

For < 220 GeV, the most visible signatures of these models are associated production 
(i/j^ip^) and pair production ijj^ijj^ ^ leading to > 3£+^t and no Z. Current LHC multilepton 
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+ analyses do not yet have the sensitivity to probe these models, due to the low and 
lepton pt when has a mass approaching its lower limit of 130 GeV, and due to the low 
production cross section for higher masses. In Section 4, however, we show the reach of such 
a search at LHC14. 

(2, -1/2) Type IV 2HDM: The particle content is the same as for Type I above, but the 
new charged particles decay almost exclusively to r leptons. Production of HA^ AH^, and 
H^H~ leads to 4r, 3t + Xt, and 2t + final states. The Z2-odd fermions decay into 
the same final states, along with extra missing energy carried away by the WIMP, x- 

The most promising channels are 
where Th is a hadronically tagged tau. The strongest constraints currently come from the 3r 
mode due to the larger production cross section. The CMS collaboration has searched for 
same-sign leptons and one hadronic tau [81] at 8 TeV and 9.2 fb""*^, and this excludes heavy 
scalar masses below 180 GeV. This is a far superior bound than one would find with other 
proposed techniques for Type IV 2HDM. In particular, we find that the significance/exclusion 
using the same-sign lepton + tau signature is approximately 4x higher for the 3t mode and 
2x higher in the 4t mode as compared with the competing analyses in [60] at fixed lumi- 
nosity! The authors of [60] find the highest significance in events with 3-4 tau tags; because 
we rely on smaller numbers of hadronic tau tags (one or two), our analysis has smaller QCD 
backgrounds and should have lower systematic uncertainties as well. In Section 4.3 below, we 
propose related searches in the same-sign dilepton + tau final states, and discuss the reach of 
such models at LHC14. Finally, we find that there are no useful constraints on the fermion 
ip in this scenario due to the limited phase space available in its decays. 

(3, 0): As in the 2IIDMs discussed above, the neutral and charged components of mix after 
electroweak symmetry breaking. With no direct coupling of (p to SM fermions, the decays of 
the scalar proceed predominantly through Higgs mixing to WW, ZZ, and hh. The scalar pro- 
duction cross section is very small, however, and such scalars are not currently constrained by 
LHC data. The production cross section of the fermion is larger, with decays tp^ W^ + x 
and tp^ h + X- However, none of the production modes is constrained by existing searches 
due to large SM top and diboson backgrounds^. 

(3, -1): The spectrum and decays are similar to the (3,0) model, with the exceptions that 
(f^ couples to ZZ and there exist doubly charged states (p , ^ . The production rate of 
the scalars is too small for observation, but the fermions decay to several highly-constrained 

^In SUSY models, the associated production of triplet fermions X2xt is strongly constrained through the 
decays X2 — ^ Zx? and xt ~^ W^Xi- 1^ such models, the neutralino decay to Z occurs through Higgsino 
mixing, whereas in the models we consider, the only allowed decay is to h, and the SUSY constraints do not 
apply. 
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final states, including ^^^^ W^Z + XX ip^^ip^ W^W^W^ + XX- The strongest 
constraints are on the W^Z + final state with leptonic decays, and 170 GeV < < 
210 GeV is excluded [80]. 

(f does not mix with the Higgs boson 

When the scalar (p does not mix with the Higgs boson it decays through its couplings 
to leptons. The model we consider involves the coupling (f e^e^. The scalar (p therefore 
appears as a same-sign dilepton resonance and is identical to a right-handed, doubly-charged 
Higgs boson . ATLAS has presented bounds on this model, and assuming 100% branch- 
ing ratio into ji~ ji~ or e~e~, if is ruled out through the entire region of interest for our 
models, m^p < 300 GeV. Decays of ^~^~X similarly excluded. However, if ^ 

decays predominantly to r leptons then the bounds are significantly weaker (m^^ > 110 GeV). 

3.4 Sterile neutrino model 

In models with a sterile state xs, the scalar messenger (f decays into a SM lepton and xs in 
Fig 6. The signature is identical to slepton production in SUSY models. The phenomenology 
depends on the mass of xs? which is a free parameter, but for concreteness we choose to 
let = for all analyses that follow. If X3 has a weak-scale mass, the constraints and 
search prospects consequently weaken. It is customary for ATLAS and CMS to quote slepton 
constraints as a function of the neutralino mass, and it is straightforward to adapt such 
constraints to the sterile neutrino model with non-zero m^g . 

Pair production of (p^ gives a dilepton+zEx signature at colliders. LEP places a bound of 
> 90 GeV, regardless of the final state lepton flavor or SU^(2) charge of (p. At the LHC, 
searches in the dielectron+zKr and dimuon+^T channels place constraints on the ip"^ mass 
provided it decays into light-flavor leptons (electrons or muons) . They exclude ip^ lighter than 
275 GeV if it has electroweak charge (2,-1/2), while (p^ is excluded between 120 — 200 GeV if 
it has charge (1,-1). There remains an open window for the SU^(2) singlet between 90-120 
GeV; while the p production rate is relatively large in this window, the events kinematically 
resemble fully leptonic WW production, a background that is two orders of magnitude larger. 
Given a mild excess already observed in WW production [85] at 7 and 8 TeV, and that the 
boosted WW fraction will rise at 14 TeV, it will be a challenge to close this gap. 

If (p^ decays to r^^-f-xs, however, there are no constraints. The leptonic branching fraction 
of taus is sufficiently small to make the dilepton-f>ET bounds irrelevant, while hadronic r- 
tagging algorithms have sufficiently large mistag rates that Z+jet events dwarf any ditau+^Er 
signal. As a result, the prospects for discovering (p^ at the LHC are very poor when it 
preferentially decays into tau, and motivates the development of high-purity tau taggers. 
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Table 3: Summary table for bounds on the scalar messenger if^ in the sterile neutrino model according 
to weak charges and flavor couplings, as well as whether ip^ can be discovered at 5cr at LHC14 with 
C < 300 fb~^. Bounds on are typically weaker in the 100-300 GeV region and are described in the 
text. 



When 99=^ is heavy (m^^ > 300 GeV), the bounds on are weak because of its small 
production cross section. We find that, in this mass range, the strongest bounds on the model 
come from constraints on the fermion tp^, whose production cross section is enhanced relative 
to (f^. The fermion ip^ decays through an off-shell (f^ into a 3-body final state X + + X3- 
As a result, bounds on ip^ are very weak when ^ m^, as much of the energy from the 
i/j decay goes into the dark matter mass. When ip^ is heavy, however, the missing energy of 
its decay is quite large due to contributions from both X3 and dark matter, and the bounds 
are correspondingly strong. If tp^ is an SU^(2) doublet decaying to electrons or muons plus 
missing energy, it is excluded from 230-415 GeV. An SU^(2) singlet is not constrained by 
current data, but bounds on the model will likely result over a similar interval at LHC14, 
while bounds on the doublet will be extended slightly down to ~ 200 GeV. The results for 
the scalar (p are summarized in Table 3. 

4. Proposals and Prospects for Future Searches 
4.1 Stable model 

As discussed in Section 3.1, the dominant way of searching for the new states ip and (p is 
via production of the charged components of the SU^(2) multiplet, and performing either a 
monojet or a disappearing-charged-track search. We mention briefly the prospects for such 
searches at LHC14. 

While some improvements in the monojet and monophoton bounds are to be expected 
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with higher cent er-of- mass energy and luminosity, the current analyses are limited by system- 
atic effects. Suggested improvements to such analyses (such as using the razor variables [86]) 
are predicted to modify the bounds by a factor of ~ 2, which is substantial but not sufficient 
for monojet and monophoton searches to constrain this model due to the small electroweak 
production rate of (^/ijj. 

Searches for disappearing charged tracks are more promising, as existing ATLAS searches 
at LHC7 already constrain some models with masses ^100 GeV [79]. An improvement in 
the bounds by a factor of two in the — '0^ mass splitting could result in < 200 being 
ruled out for an SU^(2) triplet. As such searches are still in their infancy, an increase in cross 
section and luminosity in the 8 TeV run, along with a better understanding of systematic 
errors, could be sufficient to yield such an improvement in the bound; however, we refrain 
from speculating as to what the expected improvement will be. On the other hand, SU^(2) 
doublets have lifetimes that are an order of magnitude shorter due to the larger mass splitting 
between lightest charged and neutral states; there is no bound achievable in the foreseeable 
future, and a great improvement in disappearing track searches would be necessary to probe 
such models. 

4.2 Odd scalar model 

The models that are unconstrained by current searches are those in which the charged fermion 
mixes with r leptons, and those where the new states are SU^ (2) singlets with correspondingly 
suppressed production cross sections. The Z2-odd scalars themselves are also unconstrained. 
They decay as Lp^ ^ + X therefore resemble sleptons decaying into leptons + a 
neutralino with mass 130 GeV. Since such scenarios are well-studied, we simply extrapolate 
from ref. [87] and quote that (2, —1/2) scalars can be discovered in the range 175-300 GeV 
at LHC14, while (1, —1) scalars and doublets decaying to taus cannot be discovered with 5cr 
significance. 

We now turn to our main study: singlet decaying to all fiavors of leptons, and doublet 
mixed with taus. We show that these scenarios can be discovered with 5cr significance at 
LHC14, with the exception of an SU^(2) singlet that mixes with r leptons. 

(2,-1/2): If the charged fermion mixes with the r lepton then the dominant final states 
for production of and/or t/^^ are t+t~+VV, where V=W, Z. Two final states give the 
distinctive signature of same-sign dileptons plus a hadronic tau: W+W~t"^t~ and W=^Zt"^t~. 
The final state also includes a dijet resonance from the hadronically decaying gauge boson. 
The dominant SM backgrounds are tt+W and WZ+jets. We include only prompt lepton 
backgrounds, while according to the recent ATLAS search for same-sign leptons [88], the 
prompt background constitutes only one-third of the total SM background for e^e^ and half 
of the total SM background for e^/i^. To account for the effects of non-prompt backgrounds 
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(such as heavy flavor decays or external photon conversion), we run three analyses where we 
multiply the background normalization by factors of 1, 2, and 3, respectively. While non- 
prompt backgrounds are not well-modeled by an overall renormalization of the background, 
they are typically softer than leptons from signal or prompt background, and a rescaling of 
the background renormalization is therefore a conservative way of characterizing their effects. 

We propose the following search, which is similar to the CMS same-sign dilepton + 
hadronic tau search [81], but with an added dijet resonance tag: 

1. Exactly two same-sign leptons with > 20 GeV and |7^| < 2.5 (electrons should have 
PT > 25 GeV) 

2. Exactly one hadronic tau with > 20 GeV and \r]\ < 2.5 

3. Xt > 40 GeV 

4. Z veto: all same-flavor lepton pairs^ should satisfy \m££ — mz\ > 10 GeV 

5. At least two jets with ^x > 20 GeV and with one pair having an invariant mass in the 
Z/W mass window (65 GeV < mjj < 105 GeV) 

We consider the cross section for a benchmark point, = 240 GeV, which is just above 
current collider bounds. The signal cross section after all cuts is 1.0 fb, while the background 
cross section from Monte Carlo is 0.18 fb at LHC14. We show in the left pane of Fig. 8 the 
dijet resonances associated with the hadronically decaying gauge boson. In the right pane 
of Fig. 8, we show the integrated luminosity needed at LHC14 for a 5a discovery of '0 as a 
function of m^, with the three curves corresponding to renormalization of the background by 
factors of 1 (bottom), 2, and 3 (top). We see that i/j can be discovered over the entire range of 
interest (m^ < 300 GeV) with 20 — 70 fb~^ of data, depending on the effects of non-prompt 
backgrounds. 

(1,-1): Because of the low production cross section and branching ratio to charged leptons 
relative to higher multiplets, SU^(2) singlet fermions ^ are not currently constrained. We 
propose a four- lepton search that allows discovery of < 260 GeV with 300 fb~^ when 
mixes with light-flavor leptons. When i/j mixes with taus, however, we find that even with 
aggressive cuts, the signal-to-background ratio is still too small for the search to be viable. 

The process we study is the production of V^"^V^~, where one fermion decays to Wu 
(branching ratio ^ 0.7), while the other decays to Z£ (branching ratio ^ 0.3). The gauge 
bosons are required to decay leptonically, and the final state is four leptons plus ^x- A 
distinctive feature of this final state is the existence of a Z£ resonance, and we exploit this 

^We perform a Z-veto on all same-flavor pairs (and not just opposite-sign same-flavor pairs) to suppress 
backgrounds from charge misidentification. 
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Figure 8: In the left pane, we show the dijet invariant mass (mjj) distribution in the doublet Z2-odd 
scalar model where ip mixes with taus. The production modes at LHC14 are pp ip^ip^ /ip^ip 
jj + £^£^ + Th + All cuts were applied except those involving mjj. The plot on the right is 

the luminosity for the 5cr discovery of tjj at LHC14 as a function of fermion mass. The three curves 
correspond to background normalizations corresponding to 1 (bottom), 2, and 3 (top) times the value 
from Monte Carlo to account for non-prompt lepton backgrounds. The plot begins at 240 GeV, which 
is the current bound. 



by cutting on the three-lepton mass, M^/^/^, where £^£^ are the leptons from the Z decay. We 
scan over and apply the following cuts: 

1. Exactly four leptons with > 20 GeV and I77I < 2.5 (electrons should have > 25 
GeV) 

2. Exactly one pair of same-flavor leptons with \mu — mz\ < 10 GeV 

3. Xt > 50 GeV 

4. At least one triplet of leptons with |M^/^/^ — m^\/m^ < 0.2, where m^/^/ reconstructs 
the Z mass 



The main SM backgrounds for such a process include diboson (ZZ) and triboson (ZWW) 
processes, where the gauge bosons decay leptonically. Some of the lepton pairs could come 
from off-shell photon final state radiation instead of a Z, and we also include this internal 
photon conversion to dileptons. We have checked our background Monte Carlo against the 8 
TeV CMS searches for four leptons with missing energy [89]. 

In Fig. 9, we show the trilepton invariant mass M£/£/£ distribution (both combinations) 
after cuts 1-3 at LHC14, as well as the luminosity for a 5a discovery. For a benchmark point 
of = 180 GeV, the signal cross section after all cuts is 0.2 fb, while the SM background is 
0.05 fb. Such a search allows for the discovery of t/j^ in the mass range up to 260 GeV with 
0(300 fb-^). 
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Figure 9: In the left pane, we show the trilepton invariant mass M^'^'i distribution at LHC14 as 

described in the text for the singlet Z2-odd scalar model, with production modes pp 

4£-\-^T' On the right is the luminosity needed for 5<j discovery at LHC14 as a function of the fermion 

mass. 

We attempted a similar analysis for models where ^ r^Z, and the final state from 
pair production is WZru. Requiring that there be one hadronic tau in the event, the 
backgrounds become substantially higher because of fake taus from diboson+jet events. After 
similar cuts, we find S/B ^ 0.02. Systematic errors on the background estimates dominate 
and prevent setting any meaningful bounds. 

Summary 

We summarize the current constraints and prospects for detection of the Z2-even fermions 
with various SM charges as a function of in Table 4. We organize the table according to 
the SM charges and final state flavors for each model, as well as the status of the signal. 

4.3 Odd fermion model 

Lp mixes with the Higgs boson 

(2, -1/2) Type I 2HDM: In a Type I 2HDM, the doublet scalar decays through mixing 
with the Higgs, while the Z2-odd '0 decays to dark matter and a (possibly off-shell) W/Z//i. 
We focus on the phenomenology of '0 decays, as they are a unique prediction of our model 
and are not present in a generic 2HDM. The charged fermion, '0^, decays exclusively to 
W(*)x5 while decays to both Z(*)x ^ind /ix, with the Z mode dominating up to — 
ruh + ^ 260 GeV, at which point both modes become equally important (in the range 
260 GeV < < 300 GeV, the branching fraction to Zx is about 60%). We discuss the 
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ijj Charge 


Flavor 


Constraints 


Reach (5cr@14TeV) 


Dominant Signal(s) 


Already excluded 






> 350 GeV 


N/A 


xP^^P^ i+£-W+Z ^ M + 

T7( 


(3,0) 




> 350 GeV 


N/A 




(3,0) 


T 


> 330 GeV 


N/A 




(3,-1) 




> 500 GeV 


N/A 


> 3£ plus (multiple 
channels) 


(3,-1) 


T 


> 400 GeV 


N/A 


i^£^ plus T/i (multiple chan- 
nels) 


Proposed search 


(1,-1) 




none 


< 260 GeV 
(£ ^ 300 fb-^) 


(no Z) 














T 


> 240 GeV 


< 300 GeV 
{£ ^ 40 fb-^) 




Very low sensitivity 


(1,-1) 


T 


none 


none 


^ W±z/Zr^ ^ 3£ + 



Table 4: Summary table for Z2-even fermions with the indicated electroweak charges and flavor 
couplings. When relevant, the mass reach indicates the highest that can be discovered in our 
window of interest (up to 300 GeV) and the luminosity required at LHC14 for a 5cr discovery. 



phenomenology of the scalars at the end of this section. 

We focus on two kinematic regions: when < 220 GeV and the gauge boson decays 
are off-shell, and > 220 GeV, when the decays are on-shell. The signature with the best 
prospects for discovery in the lower-mass region is the production of i/j^il;^ M+ This 
is similar to neutralino production in SUSY models, where the heavier neutralino decays into 
an off-shell Z boson and the LSP. The dominant SM background is diboson production (ZZ), 
and we also include triboson (WWZ), both with fully leptonic decays. The backgrounds are 
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Figure 10: We consider the Type I 2HDM in the ZZ + final state. In the left pane, we show the 

distribution at LHC14 after all cuts (except the cut). The production mode is pp ip^ip^ 
4£ + On the right is the luminosity needed for a 5<j discovery at 14 TeV is given as a function of 
m^. This search mode targets < 220 GeV. 



similar to those used in the charge (1,-1) Z2-even fermion search in Section 4.2 above. For 
< 220 GeV, the backgrounds are efficiently suppressed with a Z-veto. We apply the 
following cuts (which are similar to current multilepton analyses): 

1. Exactly four leptons with pT > 10 GeV and I77I < 2.5 (leading lepton has pT > 20 GeV) 

2. No pair of same-flavor leptons with \mu — mz\ < 15 GeV 

3. To remove leptons from photon conversions, we veto events where three leptons satisfy 
\Mm - mz| < 15 GeV 

4. Xt > 50 GeV 

The pt cuts are relaxed relative to our earlier M analysis because backgrounds are very small, 
and discovery potential is consequently optimized by maximizing the signal efficiency. We 
present the 0t distribution after cuts 1-3 and the luminosity needed for 5a discovery at 
LHC14 in Fig. 10. For a benchmark point of = 200 GeV, the signal cross section after 
cuts is 0.11 fb, while the SM background is 0.017 fb. Discovery can begin for ^ 180 
GeV with an integrated luminosity of ^ 80 fb~^. Discovery of ip below this mass is more 
challenging because most of the final state energy is taken by the dark matter mass and there 
is little available phase space to pass the kinematic cuts. At ^ 150 GeV, there is a feature 
in the curve of luminosity required for discovery in Fig. 10. We have confirmed that this is the 
result of the convolution of the falling production cross section, a rising efficiency of passing 
some kinematic cuts (such as lepton p^ and Xt), and a falling efficiency of the Z-veto at 
higher masses. 
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Figure 11: We show the luminosity needed for a discovery in ip^ip^ W^Z + at 14 TeV in 
the Type 1 2HDM. This search mode targets > 220 GeV. 

For > 220 GeV, the Z in ip^ decay goes on-shell, and ZZ backgrounds dominate over 
the signal A more effective search strategy considers the associated production mode 
W^Z + XX' We now require that two leptons reconstruct the Z, and the missing energy from 
XX ill the final state allows for signal discrimination through a cut on the transverse mass 
(mx = ^J^^--^) of the missing energy and the lepton not associated with the Z. We 
propose the following cuts: 

1. Exactly three leptons with px > 10 GeV and I77I < 2.5 (leading lepton has px > 20 GeV) 

2. No hadronic tau with px > 20 GeV 

3. One pair of opposite-sign, same-flavor leptons with |m^^ — mz| < 15 GeV 

4. Xt > 50 GeV 

5. mx > 120 GeV 

The dominant background is W^Z. We present the luminosity needed for 5cr discovery at 
LHC14 in Fig. 11. For a benchmark point of — 250 GeV, the signal cross section after cuts 
is 0.8 fb and the background is 3.4 fb. Discovery is possible for 230 GeV < < 265 GeV; 
for lower masses, there is insufficient energy to pass the stringent mx cut, while at higher 
masses, signal/background discrimination is, in principle, possible but the rate is too small. 

The scalars i7, A, and are also abundantly produced, although their discovery 
prospects are dim, as we now discuss. These scalar decay to b quarks, gauge bosons, and 
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the Higgs boson /i, depending on their respective masses. In general, these signatures are 
difficult to disentangle from QCD backgrounds, and we touch only briefly on the possibilities 
for detection. Below two times the mass of the bosons, the most distinctive signatures 
are HA Ah with a leading-order (LO) cross section of 235 fb for ttih — tua = 100 GeV. 
This is still about three orders of magnitude below the bb background, and observing the 
double-resonance structure of the signal events is still challenging to distinguish from the 
massive background. 

At higher masses uih = tua ^ 200 GeV, other decay modes are available to the heavy 
Higgs, H ^WW, ZZ while A still decays to bb. The final state WWbb looks exactly like tt, 
but with a cross section four orders of magnitude smaller, and this also is effectively invisible. 
Similarly, the signal ZZbb can give distinctive M + 2b signatures, but suffers significantly from 
branching fractions (leptonic and H which is ^ 20%). Indeed, we find that this is 

about an order of magnitude smaller than the background for 4 leptons from Z bosons and 
> lb jet. Finally, at masses above 250 GeV, we get the final state HA hhhZ. This can 
give spectacular 2i + 6b signatures, but the cross section of such processes is minute (~ 6 ab), 
and such a search is not feasible except at perhaps the very highest integrated luminosities 
of LHC14. 

The prospects are not any better for H^ final states. For uih < rrit^ H^ decays to ru, 
and signal events are swamped by W backgrounds. Above m^, H^ decays to tb^ giving rise 
to 2t + 26 and tbVY final states. Once again, signal discrimination from the enormous top 
backgrounds is most likely impossible without relying on extremely subtle features of the 
signal kinematics. At this point, it does not appear that the backgrounds are sufficiently 
well-understood for such channels to be viable. 

(2, -1/2) Type IV 2HDM: In the Type IV model, the charged Higgs and neutral Higgs 
decay predominantly to taus, which provide tau-rich final states. Because of the large con- 
tamination of VV+jets backgrounds to final states with leptons and taus, we concentrate on 
the relatively clean signature where one tau decays leptonically and the other hadronically. 
This leads to same-sign dilepton + hadronic tau signatures similar to those studied in Sec- 
tion 4.2. The most-constrained process is the associated production of heavy charged/neutral 
Higgs bosons, with the production of H^ + A/H ^ 3t + u. We find that the most powerful 
discriminant is the invariant mass between each of the leptons and the tau, m^^-, which is 
much larger than in the SM background due to the higher mass scales present in the signal 
production process. The SM backgrounds are ZZ and WZ. We use the following cuts: 

1. Exactly two same-sign leptons with > 20 GeV and |?7| < 2.5 (electrons have > 25 
GeV) 

2. Exactly one hadronic tau satisfying > 20 GeV and \r]\ < 2.5 
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Figure 12: In the left pane, we show the ruir distribution at LHC14 after ah cuts (except the ruir 
cut) in the Type IV 2HDM. The production modes are pp H^H/H^A £^£^ + Th + Xt- The 
mass of the heavy Higgses is set to be 200 GeV, and they decay 100% to taus. On the right is the 
luminosity needed for a 5<j discovery at 14TeV is given for different decay branching ratios to taus, 
along with three different prefactors renormalizing the Monte Carlo background: 1 (bottom), 2, and 
3 (top). 

3. Xt > 50 GeV 

4. The invariant mass of one lepton and the tau must satisfy m^^ > 200 GeV for each 
possible combination 

For a benchmark point of uih = TaA = — 200 GeV, we show in Fig. 12 the invariant 

lepton-tau mass m£r distribution after cuts 1-3, and the luminosity needed for 5a discovery 
as a function of heavy Higgs branching ratio to taus. After all cuts, the signal cross section 
for mif = 200 GeV is 2.1 fb, while the SM background is 0.06 fb at LHC14. We see that this 
search is powerful and probes 2HDM beyond Type IV: it allows for discovery of heavy scalars 
with branching ratios to taus as low as 0.2 in the LHC late running. 

There are also final states HA 4t, and the above search can be applied to this scenario 
with the following modifications: 

1. Exactly two hadronic taus^ 

2. m£r > 100 GeV for all combinations 

A more relaxed m£r cut is chosen because almost all backgrounds are eliminated by the 
second tau tag, and using a m^^ > 200 GeV cut has very poor statistics in our Monte Carlo 

^To be conservative, we multiply the signal efficiency by 50% when using > 1 tau tag to account for possible 
reductions in r-tagging efficiencies, mis- modeling of tau tagging in PGS[90], and/or enhancements in systematic 
uncertainties associated with having multiple hadronic r leptons in the final state. 
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Figure 13: We show the luminosity needed for a 5cr discovery in 3r (left pane) and 4r (right pane) 
searches at 14 TeV as a function of the heavy Higgs mass in the Type IV 2HDM with final states 
i^i^ + (1 - 2)rh +Xt. In this plot, we assume that H^,H,A decay to taus with 100% BR. The three 
curves have the same meaning as those depicted in the right pane of Fig. 12. 



study. Background would be even further reduced, but there would be significant systematic 
uncertainties associated with its estimation, and since the analysis is predominantly limited 
by the small signal statistics anyway, we consider m^^ > 100 GeV. 

In Fig. 13, we plot the luminosity needed for a 5a discovery of heavy scalars as a function 
of the scalar mass assuming 100% branching ratio to taus for the 3t (left) and 4t (right) 
final states. Although the 3t search is more powerful, the 3t and 4t analyses can serve as 
independent cross-checks of one another, and in the 3t channel with 20 — 35 fb~^ at LHC14, 
the scalars can be discovered over our entire region of interest. 

The prospects for the discovery of the Z2-odd '0 are poorer. The available phase space 
from il; decay is limited by the massive dark matter in the final state, and the neutrinos from r 
decay carry off some of the remaining kinetic energy. As a result, we find that the efficiencies 
of kinematic cuts are lower for and, in spite of the larger production cross section compared 
to H/A^ the above 3t and 4t analyses do not allow for 5a discovery of ^ in the Type IV 
2HDM with 300 fb"^ 



(3,0): The scalars in the (3,0) model decay through Higgs mixing with similar branching 
fractions as the Type I 2HDM, and are therefore challenging to discover at the LHC. 

The '0 production modes are '0^'^^ hh ^ XX? '0^'^^ W^/i + xXi '0^'^^ 
W^W^ + XX- All of these final states suffer from either large backgrounds, or are clean 
signatures but have very tiny rates due to small branching fractions (ex. hh ZZZZ ^> M). 
The most promising signature we find is W^/i + with h 77; this is still suppressed by 
the small /i ^ 77 and leptonic W branching fractions, but is one of the cleanest Higgs decay 
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Figure 14: We show the luminosity needed for a Sa discovery in ip^ilP 4£+^t searches at 14 TeV 
for the Z2-odd fermion with charge (3, —1). The shaded region is already excluded. This search mode 
targets < 220 GeV. 

modes to study. The resulting final state is ^ + 77+^X5 where the photons reconstruct ruh if 
h is on-shell in decay (i.e. > 250 GeV). This final state is similar to the conventional 
Higgs-strahlung process W^/i, although with more invisible particles, and the two processes 
have comparable cross sections at 14 TeV (^1 — 2 pb, depending on m^). 

Based on extrapolations from current W^h ^ + z/^ + 77 searches [91], it is possible 
that the '0 could be discovered over QCD and Higgs-strahlung backgrounds, particularly with 
cuts that are more targeted to the kinematics of tp^ decay, but this would likely require a high 
integrated luminosity. A precise understanding of the backgrounds, systematic uncertainties, 
and the consequent optimization of a '0 search, are beyond the scope of this paper, but the 
prospects should be studied in more detail by the experimental collaborations. 

(3,-1): The spectrum and decays in this scenario are similar to the Type I 2HDM. With 
^ 100 fb~^ of integrated luminosity at LHC14, t(; can be discovered at 5a for 100-300 GeV 
masses. For < 220 GeV, decays to off-shell V + multilepton searches with a Z- 

veto are effective. Considering in particular the production mode pp ip^ip^ Z(*^Z(*^+^T5 
we apply the M + analysis from the Type I 2HDM to with charge (3,-1). The results 
are shown in Fig. 14. Discovery of '0 can begin at ^ 15 fb~^ at LHC14, and the entire range 
up to — 220 GeV can be discovered at 5a with 50 fb~^, except in a small window where 
ip and X ^1*6 degenerate. 

For > 220 GeV, the gauge bosons from ip decay are on-shell, and there are two 
dominant production modes: ip^^lj^ W^Z + XX^ in the Type I 2HDM, and '0^^^^ 
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Figure 15: We show the luminosity needed for a 5<j discovery in ip^ip^ 3i -\- (left pane) and 
^±±^T ^ ^±^±^2j+^T (right pane) searches at 14 TeV for the Z2-odd fermion with charge (3, — 1). 
These search modes target > 220 GeV. 



W^W^W^ + XX' For the W^Z + 0t final state, the trilepton analysis proposed for the 
Type I 2HDM provides the best prospect for discovery, while the production of the doubly 
charged fermion leads to an i^i^ + 2j + Et final state. We restrict ourselves to same- 
sign dimuon signatures, since muons have small non-prompt backgrounds, and we consider 
diboson, W^W^, and ttW^ backgrounds. The cuts we propose are: 



1. Exactly two same-sign muons each with > 20 GeV and a combined invariant mass 
m^^ > 12 GeV 

2. No additional reconstructed leptons, hadronic taus with > 20 GeV, or 6-jets with 
Pt > 20 GeV (assuming a 70% tagging efficiency and 1% mistagging rate) 

3. Same-sign leptons do not reconstruct a Z (|m^^ — mz\ > 15 GeV) 

4. At least two jets with px > 20 GeV and \mjj — mw| < 15 GeV 

We show in Fig. 15 the luminosity needed for 5a discovery as a function of m^, considering 
the trilepton search in the left pane and the same-sign muon search in the right pane. We see 
that both searches can discover ^ between 220 and 300 GeV. A discovery in the ^ 230 GeV 
range is possible with 50 fb~^ with the trilepton search, while a discovery in the ^ 300 GeV 
range is possible with 20 fb~^ with the same-sign dimuon search. Both discovery modes 
complement one another and could be used to confirm any observed excess. 



(p does not mix with the Higgs boson 



One model we considered in Section 3.3 included a doubly charged scalar that decayed 
into same-sign leptons through a (f e^e^ operator. Decays of (f to electrons and muons are 
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already ruled out unless the branching ratio is < 1%. This leaves the decay of Lp into 
T~T~ . Pair production of gives a 4t final state, and this scenario can be constrained 
using the same 4r search we applied to the Type IV 2HDM above. Because the details 
are similar, we do not repeat the 4t analysis in full, but simply re-scale the signal cross 
section to match the Drell-Yan production of (p . We find that the ratio of cross sections, 
a{pp (p^^(p )/(j{pp HA)^ is approximately 1.5 at — = 100 GeV and 1.9 at 
— = 300 GeV . Applying this to the results of Fig. 13, we find that new discoveries can 
begin with ~ 8 fb~^ of luminosity, and this scenario can be completely excluded in the range 
m^— = 100 - 300 GeV at LHC14 with ^ 7b - 125 fb'^ 

Summary 

We summarize the current constraints and prospects for detection of the Z2-even scalars 
and their partner fermions in Table 5. We organize the table according to the SM charges 
and final state flavors for each model, as well as the status of the signal. 



5. Discussion and Conclusions 

We have examined the bounds and search prospects from the LHC on a variety of models 
exhibiting light (~ 100 GeV - 300 GeV) charged states. These models are directly motivated 
by the recent observation of a 130 GeV gamma-ray line in the Fermi data, although the 
resulting collider signatures apply more generally to searches for (9(100 GeV) electroweak 
states. We have classified the results according to the Z2 charge of each particle and their 
electroweak gauge charges. The results were shown in Fig. 1. 

Generally, electroweak doublets (or higher electroweak multiplets) decaying to light-flavor 
leptons with gauge bosons are completely excluded by current searches. Electroweak singlets 
and final states with multiple tau leptons or gauge bosons have weaker constraints, but 
most of these scenarios can be probed at LHC 14 with up to 300 fb~^ luminosity. We have 
proposed searches for LHC 14 that are modifications of existing multilepton analyses and that 
can discover the new charged fermion, and charged scalar, (p. In particular, same-sign 
lepton + hadronic tau searches are identified as excellent probes of 3t and 4t final states, 
and can improve the signal significance of such models by factors of 2 — 4 over other proposed 
analyses at fixed luminosity. Finally, we have found that there are some models which require 
more exotic searches (such as the disappearing charged track signature of the stable model), 
or whose signatures are completely buried in SM backgrounds (due to small production cross 
section or tt+^t final states); such models are very challenging to see in hadronic machines, 
and a lepton collider such as the International Linear Collider (ILC) may be needed to probe 
them. Overall, many models of electroweak physics generating a strong gamma ray line 
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Charge 


Flavor 


Constraints 


Reach (5cr@14TeV) 


Dominant Signal(s) 


Already excluded 


(1,-2) 




> 300 GeV 


N/A 


(p—(p++ 4£ (no Z) 


Proposed search 




N/A 


none 


< 265 GeV 


4£+Xt 


Type I 






{C = 300 fb-i) 


^±^0 ^ + XX ^ 
3£+Xt 




N/A 


mH> 180 GeV 


m,H < 300 GeV 


+ H/A -^St + u^ 
^^^^ + Th +Xt 


Type IV 






(£ « 20 - 35 fb-^) 




(3,-1) 


N/A 


^ (170, 210) GeV 


< 300 GeV 
{C « 100 fb"^) 


V'°V'° ^ z(*)z(*)+xx ^ 

4£+Xt 
3£+Xt 

■0=^=^^"^ -> 3W + XX -> 
Ai±Ai^ + 2j+XT 


(1,-2) 


T 


> 110 GeV 


< 300 GeV 
{£ ~ 100 fb-i) 


If ip++ -> 4r -> 
+ 2rh +Xt 


Very low sensitivity 


(3,0) 


N/A 


none 


none 


^±00 _^ + XX ^ 
+ 77+Xt 



Table 5: Summary table for Z2-odd fermions and partner scalars with the indicated electroweak 
charges and flavor couphngs. When relevant, the mass reach indicates the highest mass (m^ and/or 
m^) that can be discovered in our window of interest (up to 300 GeV), and the luminosity required 
at LHC14 for a 5cr discovery over the possible range. 



through dipole and Rayleigh operators are either excluded, or can be probed with moderate 
luminosity at LHC14. 
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A. Monte Carlo 

For the collider studies performed in our paper, we used MadGraph 5 [92] to simulate parton- 
level processes for signal and background, Pythia 6 [93] for showering and hadronization, and 
PGS 4 [90] for detector simulation (LHC card, anti-fcx, R = 0.5 algorithm for jets). For dibo- 
son processes and other analyses sensitive to soft radiation, we used MLM matching involving 
samples with up to two extra partons in the matrix element. Spectra and model files for the 
two Higgs doublet scenarios were calculated using 2HDMC [74] . We checked the signal efficiency 
of PGS T-tagging using rv and Z' tt events for mz'^w = 100, 200, 300 GeV; the 

efficiencies were in the range 40 — 50%, which is comparable to the benchmark values used 
by ATLAS [94]. We also measured the mistag rate on a dijet sample and found that it is 
^ 4 — 6% depending on again comparable to the ATLAS values. 

We use leading-order cross sections for signal and background for consistency, as next- 
to-leading-order corrections do not exist for all of the channels we study. The dominant 
backgrounds are electroweak, as are the signal processes, and S/B is expected to be compa- 
rable in both scenarios. 

B. The branching ratios for Z2-even fermions 

We provide details of the branching ratio calculations for Z2-even fermions with charge (1,-1); 
the other calculations proceed analogously, and we show only the results. 

The decay of Z2-even fermions with charge (1,-1) 

The electroweak Lagrangian for the SM lepton fields (£, e^) and messenger fields (V^, ijj'^) 

is: 



where H is the SM SU^(2) Higgs doublet, and we have defined sw = sin^w? etc. We 
diagonalize the mass matrix and find a left-handed charged fermion mixing angle of tan = 
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while the right-handed mixing angle is tan6^R 



-A2Ae . 

2 m: 



^; the left-handed mixing is 



largest. The mass of the charged fermion is m^± = y + (^2^;) ^/ 2 to leading-order. The 
decay widths of the fermion components are: 



r ^w- + u) 



^sin^^LCOS^^L^ 



1287r 

sm 6'l 



m 



rrir 



1 + 



w 



2 

w 



1 + 



w 



2m| 



(A2 cos cos 



^2 



m 



h 2 



647r 



, (B.2) 
(B.3) 
(B.4) 



We plot the branching ratios in the left plot of Fig. 16. 
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Figure 16: The decay branching ratios of charged and neutral fermions with electroweak charges 
(1,-1) and (2,-1/2). 



The decay of Z2-even fermions with charge (2, —1/2) 



The branching ratios are shown in the right plot of Fig. 16. 



The decay of Z2-even fermions with charge (3, 0) 



The triplet vector fermion has three components: ('0^, ) . We show the branching 

ratios in the left plot of Fig. 17. 
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(3,0) decay BR of 



(3,-1) decay BR of ij/^ 



1.0 
0.8 
0.6 
0.4 
0.2 
0.0 



1.0 
0.8 
0.6 
0.4 
0.2 
0.0 



->W +e+ 



(^"->h+v 



->Z+v 



150 200 250 

(GeV) 
(3,0) decay BR of ij/- 



i/f ->W +v 



150 



200 250 

(GeV) 



1.0 
0.8 
0.6 
0.4 
0.2 
0.0 



^o_>W-+e+ 



->Z+v 



(A"->h+v 



300 



150 200 250 

(GeV) 

(3,-1) decay BR of i/f' (i/f—) 



l/C_->Z+e _ 



300 



1.2 
1.0 
0.8 
0.6 
0.4 
0.2 
0.0 



150 



200 250 

(GeV) 



300 



i/f ->W +e 




: i/f ->w +v 




i//~->Z+e~ 


i/f ->h+e 



300 



Figure 17: 
(3,-1). 



The decay branching ratio of charged and neutral fermions with SM charge (3, 0) and 



The triplet vector fermion has three components: ('0^, ip ) . We show the branch- 

ing ratios in the right plot of Fig. 17. 



C. The branching ratios for Z2-odd fermions in 2HDMs 

In this appendix, we present and discuss the decay modes of the Z2-odd fermions with elec- 
troweak charge (2, —1/2) (Type I 2HDM). First, we note that, since the scalar in this model 
mixes with the Higgs boson and gets a VEV, the neutral fermion in turn is mixed with DM 
through the Yukawa interaction Xiptpx- The DM doublet fraction should be smaller than 
^ 0.2 — 0.45 to avoid gamma ray continuum bounds on annihilation to WW. In addition, due 
to the fermion mixing induced by the scalar messenger's VEV, the neutral fermion component 
can decay to DM + Z, while the charged fermion can decay to DM + W^. There are also 
decay channels to DM and an off-shell (f coming from the same Yukawa interaction X(ftl^X- 

The most direct decay mode is a two-body decay into DM and a gauge/Higgs boson, which 
is allowed when > 200 GeV. In this mass range, decays exclusively into W^ + x, while 
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Figure 18: In the left pane (right pane) we depict the three-body decay branching ratio of the neutral 
fermions ip^ (charged fermions ip^) with SM charge (2,-1/2) and odd under Z2. 



'0 ^ Z + X when < rrih + m^; otherwise, the Goldstone equivalence theorem predicts 
that tp^ has approximately equal branching fraction into h -\- x Z + x- 

For lower masses, ip has three-body decays, for example into lepton final states through 
off-shell gauge bosons, as we now discuss. For the neutral fermion, the decays are dominated 
by Z exchange and are shown in the left pane of Fig. 18. For the charged fermion, mediates 
similar three-body decays and the relevant branching-ratios are presented in the right pane of 
Fig. 18. The difference between rtvx ^^X the decays involving 1^* generation fermions 
is due to the contribution from scalar exchange, which is proportional to the different Higgs 
Yukawa couplings. 
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